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CORRELATION SURFACES OF TWO OR MORE INDICES WHEN THE 
COMPONENTS OF THE INDICES ARE NORMALLY DISTRIBUTED 


By GreorceE A. BAKER 


Indices are widely used in statistical analyses." In many cases incorrect 
conclusions are drawn because indices are not uncorrelated or independent even 
though all of the component variables are independent. In a previous paper’ 
the distribution of an index both of whose components follow the normal law was 
given exactly i.e. without approximation. The purpose of the present paper is 
to give the simultaneous distribution of two or more indices when each of the 
components follow the normal law. The case for two indices will be discussed 
in detail and the extension to more indices will be indicated. 

Let 21, 22, and 23, be correlated variables each being normally distributed 
about their respective means m, , m2, m3, With standard deviations o; , o2, 03, 
and let the correlations between the variables in pairs be represented by rz, 
113, T23. Then the simultaneous distribution of these three variables will be 


9 


2 


1 | — m) + Ros(a2 — me)” , Rs3(x3 — ms)” 
ens i Seen 


1 
Qr)iRioiows  —2R . 
4 Ry Ht — M2 = ™) , op, Gr — Ma) (s — ms) 
0102 0103 
4 2Ros (xz — ms)(x3 — ms) | da, dx2.dz3 


0203 
| 1 T12 «113 | 
1 123 
1923 1 | 


and R;; are the respective second order minors of FR. 


1Rietz, H. L. ‘On the Frequency Distribution of Certain Ratios,’’ Annals of Mathe- 
matical Statistics, Vol. VII, No. 3, Sept. 1936, pp. 145-153. 

2 Baker, G. A., ‘Distribution of the Means Divided by the Standard Deviations of 
Samples From Non-homogeneous Populations,’’ Annals of Mathematical Statistics, Feb. 
1932, pp. 3-5. ° 
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If we make the transformation 


Ly 
a=, Yi = 2123 
Xs 





Ze fio 
2=>—,) Lo = £223 
3 








23 = %3, T3 = 2 





dx, dx: dx3 = 23 dz; dze dzs 





which is certainly valid if x; , zz , 23, are all positive, then (1) becomes 


5 | eae, i mx)" Roo(2223 a m2)" 
-5 R 





2 + = 


G1 02 


1 
Gr)! Rea 
(2) 4 Bales —ms)* 5 op, (res — mdGres — me) 4 opp, (exze — ma) (@a — me) 


o3 0102 0193 














+ 2. ——— (e223 — ma) ———_ 23 dz, dze des. 
0203 

If x; , Z2, 2X3 are all positive the corresponding distribution of z, and 2 can be 
obtained by integrating (2) between the limits 0 and © with respect to z. 
If 21, Z2, 23 are all negative z, and z, are again both positive so that in order to 
get the total distribution for z,; and ze it is necessary to add to the integral of (2) 
between the limits 0 and « with respect to z; the similar integral of (2) with z; 

replaced by —z;. The result is 

2s 3e b 

(3) , — R' _ a oe dz + RO a 
(2r)'R? o10203| ~/2 a a /2 


where 


gw hg 4 gy Se + Oe, 4 Ate, 4 Aha 


a} a a3 0102 0193 0203 























Ru Re» Rss Ris Ry Ris 
b= — M2, + —> Meee + —F Ms + — 2M + — MB + — MX, 
C1 oa; 03 0102 0102 0103 









+ a m, + LJ M3Z2 + Be sn, 


0103 0203 0203 















. = Hn m? + Be 22 ot 4 Het 2 2h me + 2h |, mg + 2 M2™Msz. 


oO; oe a; 0102 0103 0203 
















The same result (3) is obtained for z,, and z. negative, z; positive and % 
negative, 2; negative and 2 positive. That is (3) is the simultaneous distribution 
of z, and z. The extension to more than 2 indices is immediate. The form of 
the distribution of the indices and the denominator variable is the same as (2) 
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except that a, b, and c, the coefficients of zj , z; and the constant term respectively 
in the exponent of e, will be different in that they will include the new indices and 
the exponent on the denominator variable will be the same as the number of 
indices involved. The distribution of the indices will again be obtained by 
integrating from 0 to © with respect to the denominator variable. 

The case when all of the variables x, , x2, 23 are independent is especially 
interesting. If rie, 713, 723 are all zero then R = Ry = Re = R3z = 1, Ry = 
Ri; = Ros = 0 and a, b, c, become a’, b’, c’, respectively. 


zi Z5 1 

, 1 3 

a a+ + — 
oj o2 03 


Mm 2 Me Ze ULE 


+s 


Oi 02 03 
2 2. 2 
dt ate +4 =o 
a 2 2 
oj 2 G3 


Under these conditions and the further condition that m, , mz , m3 are large with 
respect to o1 , o2 , o3 respectively so that the integral term of (3) maybe neglected 
(3) becomes 


aa , 











, +; 8 4% 
m2 m2 m2 22 22 1 
—4 enipeadioun siege lime Mm 21 4 Me Ze m3 2 
a4 4 qd o a a oe 
é e oi 02 03 
terre: Spergemnnaee ae 
zi z 1\' zi Ze 1 
1 2 1 2 
27010203 (2 + 3 + ) (3 > a) + ) 
Oo} d2 03 O1 d2 C3 


It is clear that z, and z are not independent in the probability sense for dis- 
tribution (4). 

The question as to the possibility of having the variables independent and the 
indices independent at the same time arises. Denote the distribution functions 
of 21, 22, X3, by Xi(ax), X2(xe), X3(x3) and of 2, 22 by Zi(2:), Ze(ze). Then, if 
xr; > 0,7 = 1, 2, 3 it is necessary that 


(5) [ X1(z321) Xo(z2322)X3(zs)25 dz3 = Z1(21)Z2(22) 


a and b being suitable limits. 
For instance, let 


‘ 1 i 
Xi (a) mm “a 2~ins3 
1 
" 1 
Xo(x2) = -35) 1 <2 <3 
Xe 
X;(z3) = 23, l<uw<2 
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Zi(z1) = 


Z2(22) = 


for value of z, and z, within a straight line sided area the corners of which are 
(3, 3), , 1), G, 1) and (1, 2). 2, and ze are not uncorrelated throughout their 
entire set of values but are for this particular set of values. Thus is appears 
that it is possible that the indices may be independent when the variables are, 
but not necessarily so. 


Indices should be used with care since it is very easy to draw invalid conclu- 
sions from the consideration of them. Usually it is better to use partial corre- 
lation analysis to remove the influence of a third factor than to calculate indices. 





THE TYPE B GRAM-CHARLIER SERIES 
By Leo A. AROIAN 


While much attention has been devoted to the Type A Gram-Charlier series 
for the graduation of frequency curves, the Type B series has been somewhat 
neglected. However the numerical examples to be presented later will show 
that the Type B series is very useful for the graduation of skew frequency 
curves. Wicksell' has demonstrated that the Gram-Charlier series may be 
developed from the same law of probability which forms the basis of the Pearson 
system of frequency curves. Rietz’ following Wicksell gives a derivation of the 
Gram-Charlier series based on the binomial (q + p)". Jordan’ gives a method 
for fitting Type B based on certain orthogonal polynomials which he ealls G. 
He uses factorial moments because of the resulting ease in finding the values 
of the constants. 

We shall consider the Type B series for a distribution of equally distanced 
ordinates at non-negative values of x. We shall find the values of the first few 
terms of the series and shall also show how the values of later coefficients may 
easily be found. We write the Type B series in the form 


(1) F(x) = eo + oAY(x) + c2A°Y (x) + c3A°Y (2) + cd*Y(x) + c5A°Y(x) + ceA°Y(z) 


where 


—m 


y(2) = — 
xr 


2 
m 


/ 
—s m = wi, the mean, 


(2) 
Ay(xz) = ¥(x) — v(x — 1) for zx =0,1,2,---s. 


Let f(x) give the ordinates of the observed distribution of relative frequencies, 
so that Zf(z) = 1. To determine the coefficients co , ¢) , C2, --- , ¢, we have, 
using the method of moments, 


Dlew(x) + Aw(ar) + cA*Y(x) + c:A*Y(x) +.....+ ced’Y(x)] = B(x) = 1. 
Salco(x) + aAy(x) + + ceA°p(x)] = Saf(x) = m. 
S2'[eo(x) + ady(x) +. + cA°Y(x)] = Za'f(x) = ue. 
z'leop(x) + + cA°Y(x)] = E2°f(x) = us. 
2"[eoy(x) + + cA'(x)] = E2'f(z) = ms. 
Da’ [eo (x) + + cA°y(x)] = Sa°f(x) = us. 
Ea"[o(x) + + cA’Y(x)] = Ex'f(x) = us. 
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Hence we must find the values of 


= n=0,1,2,3--- 
4 it AP ee 
@) 2PvVve, 5 o1.2.3--. 
defining A°y(z) = y(x). We assume that we are dealing with distributions 


in which s is large, and that the error involved in substituting >> x"A?y(zx) for 
z=0 


8 
> 2x"A’y(x) is negligible. To find these summations in a straightforward 
z=0 
manner would involve too much labor, so we shall briefly discuss some properties 
—m x 
‘ . em ; : 
of the generating function, ¥(x) = ——, the Poisson exponential, very useful 
x 


in the graduation of frequency distributions of rare events. The first eight 
moments about the origin are: 


wo = 1 = Ly(z), “1 = m = Lazy (z2), uw = m+ m = ra'y(z) 

ws = m+ 38m? + m = Z2'Y(z) 

wg = m + Tm? + 6m'* + m = Zer‘y(z) 

(5) ws = m+ 15m + 25m? + 10m + m = Ir°p(z) 

us = m+ 31m + 90m? + 65m* + 15m? + m® = Sa°y(z) 

uz = m + 63m + 301m* + 350 m* + 140m? + 21m*> + m’ = Za'V(z) 
us = m + 127m" + 966m" + 1701m* + 1050m° + 256m” + 28m’ + m® 





= 22°y(z) 
These may be found by the formula given by Jordan,’ 
» , dus 
(6) Moy = m( =o = 
Proof: avis) = a) y(x). 
dm m 


We multiply by x” and sum, giving (6). This result may readily be proved also 
by means of recursion formulas without differentiation. Now we must find the 
values of 


20 ee Oe ws 
Lae Ww y OR Oe 
Eee Te 
We do this by proving 
— Nas d a RNAs 
(7) Daa y(2) = — De 2"A'y@). 
z=0 m z=0 
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Now 
8) MD = yx — 1) — vl) = —Av(e). 
Hence 


Os d s 
+ (5)ve noe + |, 


since A’y(x) = ¥(x) — ()vc —1)+ (5) ve —2)+---+(-D’y@ — 8). 
Then by (8) 


Have) =[ve—-v@ —({)@—9 + (f)ve- 


+ (5) vc —3)-—- (3) vc ee 
-(-)'ve-9)|. 


@ gave =o +(°t")we-0- (3 ')ve — D+... 
— (-)"'¥@—s— 1). 
-|v@) ~ Pt are — 1) +(°3 ula —2)+--- 


+ (-1)¥@-—s—- » | 


= —A"**y(z). 


We multiply (9) by x”, sum with respect to x, giving (7). 
Thus by use of (7) and (5) we get: 


LA’y(z) = 0, p=1,2,3,--- 


ssn) = = 

D2 Ay(x) = — o re’ y(x) = ny + m’) = —2m — 1. 
(10) dm dm 

>r'Ay(r) = —3m — 6m — 1. 

D2‘Ay(x) = —4m*® — 18m? — 14m — 1. 


ra°Ay(xz) = —5m* — 40m*® — 75m — 30m — 1. 
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(10) 
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ra°Ay(x) = —6m*> — 75m* — 260m* — 270m" — 62m — 1. 
SrA°y (x) = 0, Ta Ayr) = 2, Va°A*y(x) = 6m + 6. 
Ya'A'y (2x) = 12m? + 36m + 14. 

rae Ay (xr) = 20m? + 120m? + 150m + 30. 

ray (x) = 30m* + 300m*> + 780 m* + 540m + 62. 


YrA*y(xz) = 0, Ta A*y(x) = 0, Sa'a’y(x) = —6. 
Seay (x) = —24m — 36, Yr’ dA*¥(x) = —60m™ — 240m — 150. 
rar A(x) = —120m* — 900m — 1560m — 540. 


SrAy(xr) = 0, Ya AV(r) = 0, Sa'a*y(x) = 24. 
ra A(x) = 120m + 240,  Sa2°d*y(x) = 0. 

raat (x) = 360m” + 1800m + 1560. 

SrA°y(r) = 0, Trd°y(x) = 0. 


Sr A°y(r) = 0, Da'A°y(x) = 0. 
rr'A*y(r) = 0, Ta'A°V(xr) = 0. 


Sa'ay(xr) = 0, Da'A°y(xr) = 0. 


Sr'A’y(x) = —120, =2x'd°y(x) = 0. 
r2'A'y(r) = —720m — 1800, Y2°A°y(x) = 720. 


Finally we substitute from (5) and (10) into (3), and for u, Wwe substitute 


kn = 


(11) 


n ; 
2. ( ) arm . Hence 
, 


r=0 
Co = 1 
cq; = 0 


Co = > (ue — mM). 


1 
Cs = — ay (us — 3u2 + 2m). 
C= a [us — Guz + po(1l — 6m) + 3m(m — 2)). 


Cs = nga — 10u4 — pu3(10m — 25) + 50pu2 (m — 1) — 4m(5m — 6)]. 
5! 


Cs = rate — 15us5 + ws(85 — 15m) + p3(130m — 225) + po(45m? — 375m 


+ 274) — 15m* + 130m” — 120m]. 


It may be asked whether criteria may be given as guides for the use of Type B. 
In general Type B may be tried if either the skewness of the distribution to be 
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‘ s 3 ° . ry 
fitted is considerable, a3 = —; > .6, or if m = we = ws approximately. The 


= 


latter condition strictly would mean that (x) alone is sufficient for a good 
graduation, if the fourth moment, us, is not used. The examples which follow 
are arranged to facilitate comparison with the Pearson system of frequency 
curves. We have an example each of Type I, II, IV, V, VI, and an example of 
the normal curve. 

Type I. Table 1. Here ag > .6 although m + po + ws. The first four 
moments, unadjusted, give an excellent fit by Type B, which is not quite as good 
as Type I. The degrees of freedom, according to Fisher,* have been taken into 
consideration here in applying the 2” test. The two classes 13, 14, were grouped 
together for the x° test. The actual numerical work is easily done on a cal- 
culating machine, although logarithms are necessary to find the value of e™ 
This example and the remaining are all taken from Elderton’ with the exception 
of Type IV which is from A. Fisher.° 

Type III. Table 2. The unadjusted moments are used. Here a3 = 2.0833 
> .6, and m = pe approximately. The fit by Type B is slightly better than that 
by Type III. We have for Type III P(x’ > 12.8) = .007, n = 3, while for Type 
B, P(x’) > 9.4 = .025 n = 3. Moreover the standard error of prediction for 
Type III is 11.2 and for Type B is 7.7. 

Type IV. Table3. The rough moments were used. Although a; = .48 < .6, 
Type B gives a fine fit since m = we = us; approximately. Here the results are 
given for Type B using 2, 3, and 4 terms of the series. This was done to show 
how the distribution changes with the addition of more terms. The superiority 
of Type B over Type IV is evident. The results for Type IV are taken from the 
class notes of Professor C. C. Craig. 

Type V. Table 4. Using the adjusted moments we have a comparison among 
Types V, A, and B. While the graduations may seem satisfactory, the 2” test 
shows that the fit is poor in each case. The order of merit is Type V, Type B, 
and then Type A. The negative frequencies which appear in Type B may be 
due to the use of the adjusted moments. If we use the rough moments, the 
negative frequencies disappear. On the whole the fit by means of the adjusted 
moments is superior. 

Type VI. Table 5. Type VI using the adjusted moments gives an excellent 
fit. Even though az; is considerable, and we = w3 approximately, four moments 
with Type B give a poor fit, and five moments, adjusted, achieve a very small 
gain. Five moments using the unadjusted moments give some improvement, 
but the —2 frequency in the first class is objectionable. 

Normal Curve. Table 6. The normal curve provides a fine fit. P(z’ > .9) = 
.96,n = 6. The first two and the last two classes were grouped together for the 
test. The fit by Type B is less probable, P(x” > 8) = .15,n = 5. Type B has 
two discrepancies, the negative frequencies, and the fact that the total fre- 
quencies (neglecting the —1) is 352. That Type B does so well is in itself 
quite amazing! 
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TABLE 1 





| Frequency computed Frequency given 














* | Actual frequency by Pearson Type I by Type B 
0 34 | 44 42.4 
1 145 137 121.3 
2 156 | 149 168.7 
3 145 | 142 156.8 
4 123 | 127 120.5 
5 103 | 108 | 94.9 
6 | 86 | 88 | 82.9 
7 | 71 69 72.2 
8 55 51 | 56.7 
9 37 | 36 | 38.0 
10 | 21 | 24 | 23.1 
11 | 13 14 12.0 
12 | 7 | 7 5.7 
13 | 3 3 2.4 
14 | 1 1 9 
m = 4.175 a3 = .712247 Type I P(x? > 4.36) = .88 
be = 7.66237 a, = 2.95214 n (number of degrees of 
us = 15.1069 Co = 1.74368 freedom) = 9 
wa = 173.326 C3 = — .078298 Type B P(x? > 9.67) = .37 
Ce = +.094592 n= 9 


F(x) = ¥(r) +1.74368 A*y(x) — .078298 A*y(x) + .094592 Aty(z). 
TABLE 2 














x | Actual frequency | "Frequency computed) — Fequeney, by 
0 44 | 59 | 48.1 
1 135 | 111 | 121.6 
2 45 | 45 58.5 
3 12 | 20 | 10.4 
4 8: | 9 3.5 
5 3 | 4 | 4.3 
6 1 | 2 | 2.9 
7 3 1 | 1.2 
m = 1.33466 a; =" = 2.0833: = ~—«.05356 

we = 1.44179 Me Cs = —.32510 

us = 3.60662 


F(x) = ¥(x) + .05356A°%y(xr) — .32510A%Y (x) 
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TABLE 3 
Number of alpha particles from a bar of polonium in intervals of 4 of one minute 





| | T | 
ype B Type B Type B 
* Frequency Type IV 2terms' | 3 terms | 4 terms 





49.0 | 58.2 
201. 199.8 
404. 386.1 
533 523.9 
521 532.1 
402 418.2 
254. 260.2 
136. 134.0 

63. 56.7 

26. 22. 

9 
3. 


57 50 | S49. 
203 | 183 201 
383 392 | 403 
525 | 544 532 
532 539 520 
408 | 417 | 402 
273 | 250 | 254 
139 =| 131 137 

45 | 61 | 64. 

27 | 26 26. 

10 | 12 | 9 

4 4 3 


oOoOoONOQorhwWN FH © 











SHOOK HOH DAAWRWE 
SCH OHOAONONRAMBDWO 





0 | 1 
1 | 0 
1 0 





m 3.87155 = .47844 
be 3.69477 = 3.506536 
3 3.39791 

uy = 47.86888 


F(z) = ¥(z) — .08839A*%(xr) — .00930A%Y(xr) + .16810A‘4Y(z). 


Type B, 4 terms P(x? > 4.50) = .72,n = 7 
Type IV P(x? > 10.8) = .15,n =7 
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TABLE 4 
Mortality Among Female Nominees 

















2 | Deaths |yoew | Tye A | 7¥Bme | SYerme | S¥onme | oteeme 
0 4 | 4 2) 14] -69 | -4 | 41 
1 is | 10 | 15 | 26.3 7.1 | 9.4 | 18.1 
2 53 | 80 | 78 | 109.7 | 100.1 | 84.6 77.4 
3 265 | 261 | 235 | 248.3 | 268.4 | 252.3 | 242.5 
4 | 4388 | 441 | 426 | 379.5 | 418.8 | 425.9 | 427.4 
5 | 525 | 480 | 521 | 432.7 | 461.0 | 484.0 | 494.1 
6 | 342 | 381 | 411 | 388.8 | 388.4 | 402.6 | 408.1 
7 | 253 | 247 225 | 285.4 | 263.5 | 259.0 | 253.9 
8 | 128 | 137 | 107 | 170.8 | 145.5 | 132.2 | 124.9 
9 | 82 | 68 | 66 84.3 68.3 | 58.6 | 54.1 
10 | 28 | 32 44 | 32.9 28.2 | 26.2 | 26.4 
11 | 12 | 14 22 | 8.6 11.0 | 13.9 | 16.4 
12 |; 8 | 6 8 | -—.01| 4.7 | 8.2 | 10.7 
3} 5] 8 2 —2.1 | 21] 43 |] 5.9 
14 | 1 | 1 0 —1.5 13 | 2.0 | 2.5 
Adjusted moments: Rough moments: 


m = 5.30435 az = .703564 m = 5.30435 
wo = 3.573345 as = 3.996196 vo = 3.65668 


v3 = +4.752437 v3 = 4.752437 
wy = 51.02659 vs = 52.85276 
Ms = 193.439125 vs = 197.39949 


Type A: f(t) = g(t) + .117261 ¢3(t) + .041508¢4(t) 
Type B: F(x) = ¥(x) — .86550A°Y(x) — .77352A%y(z) 
+ .02814A4Y(r) + .57459A5Y(z) 
Using uncorrected moments 
Type B: F(z) = (x) — .82384A°(x) — .73185A%*y(z) 
‘ + .03192A4)(x) + .94033A5y(z) 
(last column above) 
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TABLE 5 
£ Frequency Type VI | 5 Sa = 
0 1 1 | —9.5 —2.0 
i 56 50 | 83.2 69.9 
2 167 168 | 141.6 143.1 
3 98 100 | 102.3 110.7 
4 34 | 36s | 41.5 40.2 
5 | 9 i 8.7 4.6 
6 | 2 2 | .05 2.0 
7 1 5 | —.4 1.0 








Corrected moments: Rough moments: 


m = 2.402174 m = 2.402174 
Me = .928835 we = 1.012169 
ws = .893096 us = .893096 
a = 4.088800 us = 4.313176 
Ms = 11.28304 
ag = .87704 
agqa= 4.2101 


Type B, adjusted moments: 

F(x) = W(x) — .73667A°y(x) — .48516A*y(x) — .06424A4Y(x) + .10365A5Y(z) 
*Type B, rough moments: 

F(x) = W(x) — .69805A2y(x) — .44654A%Y(2) — .06587A4Y(x) + .15165A°Y (x) 


* This is used in last column of above. 


There is a slight error here, which however will 
not affect the results materially. 


The third decimal place may be slightly wrong. 
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TABLE 6 
Normal curve 


8 


Frequency Normal curve 
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Moments corrected: 


m = 5.393443 

He = 2.769635 

ws = .029805, uy = 22.40663 
az; = .0064 

a, = 2.920997 


Type B: F(x) = (x) — 1.3119A%(x) — .4179A%Y(x) + 2.1625A4Y(z) 
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A TEST OF A SAMPLE VARIANCE BASED ON BOTH TAIL ENDS OF 
THE DISTRIBUTION 


By Joun W. FeErtTiG 


WITH THE ASSISTANCE OF ELIzABETH A. PRoEHL’ 


(1) Introduction 


In testing the hypothesis, say Ho, that an observed sample E of size N has 
been drawn from a normal population for which the standard deviation, o, has a 
particular value, oo , one may form the ratio 


, Nad 
v= S (x; — m)*/o5 = —= 
i Go 


if the population mean m be known, or 


N 2 

i Ns 

v = S (x; —2)*/o5 = — 
i=] Oo 


where < is the sample mean, if the population mean be unknown. The proba- 
bility of obtaining a larger (or smaller) value of v or v’ than that observed may 


readily be obtained from the appropriate tail area of the x’ distribution with 
n = N or n = (N — 1) degrees of freedom respectively. The alternative 
hypotheses to Ho concerning the normal populations from which the sample 
may have been drawn assign different values to o and form a set of hypotheses, 
2. The members of 2 may be classed according to whether they specify 
¢>o0),0re < oo. The practice of regarding only one tail of the distribution, 
the upper or lower depending on whether v > N or v < N, is tantamount to 
accepting as admissible alternatives to Hy only one of the classes of 2. 

The alternatives may sometimes be limited to one class or the other through 
some a priori knowledge, or the problem may be such that only one of the classes 
is relevant. However, since this is not generally the case, some method of 
considering all of the alternatives is needed. When testing hypotheses con- 
cerning the mean of the sampled population, the problem is quite simple, since 
the distribution of means is symmetrical. Thus, the “corresponding” value to 
any positive deviation, ( — m), is the negative deviation of the same magnitude. 
Merely doubling the tail area pertaining to either of the deviations will serve to 
take account of both classes of alternatives, i.e., those in which m > mp and 
those in which m < mp. The problem is more difficult in the case of v or v’, 


1From the Memorial Foundation for Neuro-Endocrine Research and the Research 
Service of the Worcester State Hospital, Worcester, Massachusetts. 
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since the distribution is not symmetrical. In addition to the value of v or 9’ 
pertaining to the observed sample we require a “corresponding” value at the 
other end of the distribution. The definition of “corresponding” which is 
accepted will determine the required value. There may be a number of such 
definitions but not all of these will be equally acceptable. The value of » 
which delimits an equal tail area specifies one of the possible definitions of 
“corresponding.”” Another definition would require that the ordinates at the 
two values of v be equal. 


The Neyman and Pearson Approach. Generalized procedures for testing 
statistical hypotheses have been elaborated in recent years by J. Neyman and 
E. 8. Pearson (1-5). These have considerable philosophical appeal and will be 
traced as a basis of solution of the immediate problem. A test of a hypothesis 
Ho consists essentially of a rule for rejecting Ho when the observed sample E 
falls within a suitable critical region w of the N-dimensioned sample space W, 
and of accepting Hy when E£ falls in (W — w). In testing any hypothesis two 
types of error may be made: 

i) Ho may be rejected when it is true; 

ii) Hy may be accepted when some alternative hypothesis, H;, is true. 
Errors of the first kind may be considered “‘equivalent”’ since, if a true hypoth- 
esis is to be rejected, it is immaterial which one is chosen. Furthermore, the 
first type of error can be controlled through our choice of the size of w, say a. 
The size cf w represents the probability of a sample EF being an element of w 
when the hypothesis Ho is true. This probability may be designated briefly as 
P{Eew|Ho}. Then 


P{E ew Hp} -|-- [ r@ | Heyday --- dey = a peaewe (IIT) 


where p(E | Hy) is the elementary probability law of the sample when Hp is 
true, i.e., 


P(E | Ho) = p(t, 22, --- tv | Ho)..........0005 (IV) 


Errors of the second type, however, are not equivalent, since their consequences 
depend on the difference of the true hypothesis from Hy). The utility of a test 
of Ho will depend largely on how it controls the second type of error. Ideally, 
the selection of a critical region should take into consideration the probabilities 
4 priori of the hypotheses composing 2. Since these probabilities are generally 
unknown, tests may be sought which are valid independently of them. 

A distinction must be made between simple hypotheses which specify com- 
pletely the elementary probability law of the sample, p(Z), and composite hy- 
potheses which specify the law subject to one or more undetermined parameters. 


(2) Simple Hypothesis Concerning Population Variance 


A test based on a critical region wy may be called independent of the probabili- 
ties 4 priori of the alternative hypotheses if it is more powerful than any other 
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equivalent test for all of the alternative hypotheses (3). An equivalent test 
is one based on a region w, of the same size, a, i.e., 


P{E €wo| Ho} = P{Eew,| Ho} ea ckewan awe eke (V) 


The power of a test based on any critical region, as w, , is the probability of its 
rejecting a hypothesis Hy) when some other hypothesis H; is true. That is, 
it is the probability of E falling in w, when H; is true. Denote this power by 
P\Eew,|H;}. The greater the power of a test, the smaller the risk of the 
second type of error. If tests as defined above exist, they minimize the proba- 
bility of the second type of error. Furthermore, the probability of the first 
type of error is no larger than a. Neyman and Pearson (2) have designated 
regions satisfying this definition as Best Critical Regions for testing H» with 
regard to the set 2. If there is no such Best Critical Region, some compromise 
region must be chosen. 

A necessary and sufficient condition for wo to be a Best Critical Region with 
regard to an alternative H; is that within w» 


P(E | Ho) < kp(E|H))........00ceeec ee ee. (VI) 


where k is some constant depending on a. If this inequality is true for any H;, 
wo will be a Best Critical Region for the set Q. 

Neyman and Pearson (2) have shown that in testing the hypothesis that 
¢ = oo, when the population mean m is known, there are two Best Critical 
regions, one pertaining to the class of alternatives for which o < oo and defined 
by v < w,, the other to the class ¢ > oo defined by v > wv. v, and v2 are values 
of v so chosen that the size of the critical region shall be a. Although there is 
no Best Critical Region for all of the alternatives, the choice of a compromise 
critical region should still depend on its control of the second source of error, 
that is, on its power for the various alternatives (4). Such a compromise 
region may be designated as a Good Critical Region. What is needed is a 
region Wy of size a defined by the inequalities v < vy, andv > vw. If » and v 
are taken as the values cutting off equal tail areas, then the power of the test 
will be less than @ for some values of o less than oo. For those values of o, Ho 
would be accepted more frequently than if it were true. Thus a first require- 
ment for a Good Critical Region is that its power should nowhere be less than a, 
the value when Hp is true. Of all such unbiassed Critical Regions of size a, 
wy should then be selected so that its power is everywherd greater than that of 
any other equivalent unbiassed region. 

Critical Regions sufficiently satisfying the above requirements can often be 
obtained by stipulating that the first derivative of the power function with 
respect to 6, the parameter under consideration, shall be zero at 6 = 6, and 
that the second shall be a maximum there. Then not only does the probability 
of the second source of error decrease as we move away from 4 , but it decreases 
most rapidly in the vicinity of 6). Critical Regions satisfying these conditions 
are called unbiassed Critical Regions of Type A, (4). Under certain assumptions 
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concerning the nature of the elementary probability law p(E | @) it can be shown 
that wo is defined by the inequalities g; < ¢; and ¢; > cz where c; and ¢ satisfy 
the conditions 


I : p(¢r) dg, eee ee aio, Rs (VII) 

f $1 p(¢r) dy PF aeGccu ea reer eas (VIII) 

where gi = OMEREIO iiveisseetes (IX) 
dé 6=6, 


and p(¢;) is the distribution function of ¢ . 
In applying these results to the testing of the hypothesis that ¢ = oo when 
the population mean is known, 


yi = (v —wN / pete ad eran uate ot Sie Saeae ais be, sue sete (X) 


Obviously p(v), the distribution of v, may be considered instead of p(¢1). u is 
defined by the inequalities v < v; and v > ve where 


I p(v) dv + I p(v)dv=mta=a............ (XI) 
0 v2 
| (v — N)p(v) dv = ve"? | = 0... (XII) 


wo so defined is also of type A;, that is, its power curve lies everywhere 
above that of any other equivalent region, vanishing in the first derivative at 
o = oo, (4). 

The use of wo as the appropriate critical region is equivalent to the use of r 
as a test criterion, where 


pr? ot cs r) ere ae ee ee ee ee ee ee a ee XIII) 


That is, a value of v yielding the same r as the observed v may be taken as the 
corresponding value. Reference to the appropriate tables and summing of the 
two tail areas gives P,, the probability of obtaining a smaller value of r when 
Hy is true. Ho may be rejected if P, is less than some previously fixed number, 
say a. If the distribution of r could be evaluated the necessity of dealing with 
two values of v would be obviated. 

The criterion r is equivalent to that deduced by the use of maximum likelihood 
ratios (6). Thus, 


N 


sia — 8 (ei—m)?/20? rT 
p(E | 0”) - (Qra*) *'*¢ eens eee Vee es (XIV) 


2 The solution is the same in terms of o?. 
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ie 2 - 2 
Maximizing p(E | o°) for fixed E and all possible o° we have 


Dmax.(E |o°) = nv"! oe s (x; — nm | a. ee (XV) 





p(E | o5) —N/2.N/2_ —}(v—N) y 

ae ere ae cows. XVI 
Pmax.(E | o*) on 
Te kin veuicw a annnncea Busine snaes (XVII) 





6 8 Lo . 12 14 6 is° 


Fic. 1. Graph of Equation xz — log, t = k log. 10 


The h* moment coefficient of \ about zero, u;,(d), is given by 


ea 


un(d) oe T(N/2) 


(Qe/NY""* (a + APO”... (XVIII) 
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For N infinite, (—2log,d) will be distributed as x” with one degree of freedom. 
For finite values of N, however, we have not been able to evaluate the dis- 
tribution of \, although the distribution of the Incomplete Beta Function serves 
as a good approximation. Approximate distributions for several values of N 
have been obtained. P),, the probability of obtaining a smaller value of \ 
than that observed, as obtained from these distributions agrees well with the 
sum of the tail areas pertaining to v; and ve yielding the same value of \ (or r). 
The construction of tables is simplified by taking (1) 


logio \ = N/2(logioe — k) 
That is, 
ee (XX) 


where « = v/N. Equation (XX) is independent of N and may be solved once 
and for all for z, given k.’ In Figure 1 is plotted the graph of equation (XX). 
For convenience, the branch of the curve giving the roots greater than unity 
has been folded back with altered scale from the minimum value of k, logwe, 
occurring at x = 1. Table I was then constructed by multiplying the two 
values of x for a given k by (N/2), referring to the Tables of the Incomplete 
Gamma Function (7) with p = (N — 2)/2, and adding the resulting two tail 
areas. The values for the odd numbers above 12 were obtained by interpolating 
between the even numbers. For N = 1, (x)? was used as a normal deviate. 
The values in Table I should be correct to four decimals. Table I is entered 
with the number of degrees of freedom, n, en which z is based. In the case of the 
simple hypothesis this is NV. , 

The following may serve as an illustration: Blood urea nitrogen determinations 
(mg./100 ce.) were made on a sample of 25 schizophrenic patients. The mean 
was found to be 15.56, the variance, 10.486. Previous investigation of blood 
urea nitrogen on a large sample of normal control subjects gave a mean of 16.03 
and a variance of 20.268, which for the purpose of the example may be considered 
as the population parameters. Then we may wish to test the hypothesis that 
the variance of the sampled population, o” , is o> = 20.268, knowing the mean 
of the sampled population to be 16.03. Calculate 


_ 8 +@—m)" _ 


a5 








.528 


Referring to Fig. 1, the value of k is about .505. Turning to Table I with 
k = .505,* n = 25, P is found to be .0457. We should thus be inclined to reject 
the hypothesis. 

For N small, the area of the tail of the distribution near zero is considerably 
larger than that at the upper end. As WN increases the distribution of v becomes 


3 Tf the solution were explicit the distribution of \ could easily be deduced from that of z. 
‘ k obtained directly from (XX) is .507, corresponding to P = .0427. 
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more and more symmetrical and the two areas approach equality. Even for 
N = 50, however, they are rather unequal, so that merely doubling the area 
pertaining to the observed v does not give a sufficiently accurate approximation. 
For N > 50 an approximation correct within several units in the third decimal 
place may be obtained by taking ~/2N(+/z — 1) as a normal deviate. This 
assumes that the standard deviation is normally distributed with variance o3/2N, 


(3) Composite Hypothesis Concerning Population Variance 


Here Ho specifies only the value of the parameter 6 = 6 , leaving undetermined 
the value of a second parameter, v. Thus, Ho consists of a subset, w, of simple 
hypotheses, each of which specifies a different value for v. Any simple hypoth- 
esis specifying different values of both parameters, 6 and », is an alternative 
to Hy. These alternatives form the set 2. The elementary probability law 
determined by Ho is p(E | Ho) = p(E | %v), while that determined by an alterna- 
tive hypothesis H; is p(E | H:) = p(E | 6;). In testing composite hypotheses 
the first requirement is to find regions “similar” to W with regard to », i.e., such 
that the chance of rejection of a true hypothesis, P{E ew | Ho}, equals a@ for all 
the values of v specified by the simple hypotheses composing Hy). A test based 
on a similar region wo) may be called independent of the probabilities 4 priori, 
if its power with respect to all the alternatives of 2 is greater than that of any 
other similar region w, of the same size, a, (3). Let 


yg. = 0 log p(E | 6v) / av lo=0, pice Pak tga ecco Nich ah ava (X XI) 


Then the equations gz = constant will describe hypersurfaces in N-dimensioned 
space, on one of which the observed EF must fall. Under certain assumptions 
pertaining to the law of elementary probability it can be shown (2) that a 
necessary and sufficient condition for w to be a similar region is that 


P\E €w(¢e) | Ho} = aP{E « W(¢e) | Ho}........... (XXIJ) 


for all values of ge , where w(¢2) and W (ge) are parts of the surface g. = constant 
common to w and W respectively. <A similar region is then built up of these 
parts w(¢2) obtaining for the various values of g.. The Best Critical Region, 
wo , for a particuiar simple alternative, H;, must then be composed of pieces, 
wWo(¢2), maximizing P{E € wo(g) | H;}. The problem is the same as for simple 
hypotheses except that we shall be working in a space W(g2) of (N — 1) dimen- 
sions. wWo(g2) is defined by the inequality 


p(E| Hi) > kg) p(E\ H)...........--0- (XXIII) 


where k(g2) is some constant depending on a. If wo(ge) is the same for all H;, 
then wu» is the Best Critical Region for testing H» with respect to ©. 

Neyman and Pearson showed (2) that in testing the composite hypothesis that 
o = oo When the population mean is unknown there are two Best Critical Regions 
corresponding to the class of alternatives ¢ < oo and ¢ > oo, defined respectively 
by the inequalities v’ < v, and v’ > vy. If the whole set of alternatives, @, is to 
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be considered some compromise region must be sought. Dealing with the case 
where similar regions exist Neyman (5) defines a Critical Region as unbiassed 
and of Type B if the first derivative of the power function, ?(E ew | H;), with 
respect to @ vanishes at 6 = 6), and if the second derivative at that point is a 
maximum. Let 


a log p(E | 6) | 
qa = —. Gree dnb en ine anand 2 
06 6=65 


Then it can be shown that the desired region will be defined by the inequalities 
gi < kilge) and g: > ke(ye) where ki(¢2) and ke(ye) are determined to satisfy 


k2(¢2) 
[ P(oige) dg, = (1 — a)p(gr)............. (XXV) 
v1 @2) 
and 
ko(¢2) oo 
i , gip(¢i¢:) dg, = (1 — a) giplgigr) dgi...... (XXVI) 
ril¢e2 —co 


where p(¢e) is the distribution function of ge, and p(¢yige) is the simultaneous 
distribution of g; and ~. 

Applying equations (XXV) and (XXVI) it follows that the appropriate 
Critical Region is defined by the inequalities v’ < v; and v’ > vy where 


vy] o 
a=ata= | pv’) dv’ + i p(v’) dy’ .......(X XVII) 
0 v 
and 
fea HF nk cneceesees (XXVIII) 


where p(v’) is the distribution function of 0’. 
The use of the unbiassed Critical Region of Type B corresponds to adopting 
as a criterion 


(N—1)/2 


u Cl Mc caincnvenwanpenans (X XIX) 


Since v’ derived from a sample of size N is distributed as v derived from a sample 
of size (N — 1), it follows that 7’ is equivalent to the r of eyuation (XIII) based 
on a sample of size (NV — 1). Therefore Table I may also be used for testing 
the hypothesis that ¢ = oo whatever be the population mean, by entering with 
the number of degrees of freedom, N — 1. 

In the example previously used, compute 


From Figure 1, k is approximately .51, corresponding to P = .0422. 
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r’ is not the same as the maximum likelihood ratio \’ (6). 


7 2 
a Pmax (E | om) _ at—N/2, pN/2,—4(0'—N) 


= ~ — py” e Z sa N” N *—* 2 yf yr’. a Xxx 
Pmax (E | o”m) ( ) 


As N becomes infinite the distribution of \’ is the same as that of the \ of (XVI). 
For N = 49, the probabilities corresponding to \’ agree with those using r’ to 
within a unit in the third decimal. 

The ’ test is biassed as may be seen in Figure 2 where we have plotted the 
power of the test based on the region w defined by v; = 3.187, v, = 22.912 for 
which a = .0436 + .0064 = .0500, on the assumption that a=. 1.0, for N = 10. 
Although the criterion is biassed it is slightly more sensitive to alternatives 
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Fic. 2. Comparison of Critical Regions for v’. Ho Specifies «5 = 1.0. N = 10. 


specifying o < 0 than is the unbiassed Critical Region of Type B defined by 
vb} = 2.953, vy = 20.305, a = .0339 + .0161 = .0500. The criterion of con- 
stant distribution, p(v’), 


N-3)/2_—4v’ yw 
ai iin bib 6 SKS oa eees (XXXI) 
. . ’ , ’ - si 
has also been considered. In this case v; = 1.903, vo = 17.391, a = .0071 + 
.0429 = .0500. This criterion is biassed for some alternatives specifying 


o < oo, but its power curve lies above that of the unbiassed region for o >a. 

Apparently the bias may be shifted at will by changing the exponent of v’. 
This may be desirable if greater weight is to be given to one class of alternatives. 
In fact decreasing the exponent of v’ to 0 produces the Best Critical Region 
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for the class of alternatives specifying o” > o, and defined by », = 0, v2. = 16.919 
for a = .0500. No region can be found giving greater power. On the other 
hand this region is insensitive to alternatives of the other class. Increasing the 
exponent indefinitely produces the Best Critical Region for the other class 
defined by ve = © and»; = 3.325 for a = .0500. 
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ON THE POLYNOMIALS RELATED TO THE DIFFERENTIAL EQUATION 
ldy_ mtart  _N 
ydx bo + bx + box? D 
By FRANK 8S. BEALE 


Introduction. In a previous issue of this Journal,’ E. H. Hildebrandt has 
established the existence of a general system of polynomials P,,(k, x) associated 
with the solutions of Pearson’s Differential Equation 


ydx D 
N and D being polynomials in z of degrees not exceeding one and two respectively 
with no factor in common. 

It was shown that the polynomials P,(k, x) = P, themselves satisfy certain 
differential equations and a recurrence relation. The classical polynomials of 
Hermite, Legendre, Laguerre, and Jacobi are special types of P,(k, x). Since 
the classical polynomials are employed rather extensively in statistical theory, 
certain of their properties are of special interest. 

It is the purpose of this paper to determine from Hildebrandt’s general equa- 
tions some new properties of P,(k, x) and to apply these properties to the 
classical polynomials. The paper consists of two parts. In part I some 
theorems are established concerning common zeros of Dand P,. In particular, 
a theorem is established to exhibit the conditions under which the zeros of P,, 
which are not zeros of D, are simple. In part II a method is outlined for the 
classical polynomials by which one can determine the number and location of 
the real zeros in the various segments into which the zeros of D divide the z axis. 
The points of inflexion and the degree of the polynomials are also considered. 

A new feature of the method employed is, we believe, its being based upon the 
use of differential equations of first order, for most part, while other investi- 
gators’ have employed differential equations of second order. As to the results 
obtained, the author believes.them to be partly new. They have points in 
common with the results of Fujiwara, Lawton and Webster. 


1 Systems of Polynomials Connected with the Charlier Expansions, etc., Annals of Math. 
Stat., Vol. II, 1931, pp. 379-439. 

2M. Fujiwara: On the zeros of Jacobi’s Polynomials, Japanese Journal of Math., Vol. 2, 
1925, pp. 1, 2. 

W. Lawton: On the zeros of Certain Polynomials Related to Jacobi and Laguerre Poly- 
nomials, Bull. Am. Math. Soc., Vol. 38, 1932, pp. 442-449. 

M.S. Webster: Thesis, Univ. of Penna. These results were kindly communicated to 
me by Dr. Webster. 


206 





ON THE POLYNOMIALS RELATED TO PZARSON’S DIFFERENTIAL EQUATION 207 


I. Theorems Concerning Common Zeros of P,(k, x) and D 


The following equations will be employed later: 
Pasalk, x) = [N + (k — n)D'|P,(k, x) + DP,(k, 2). 


=~» "| P,(k, 2). 


Pislk, 2) = (n+ p | os 


Prailk, x) = (N + (k — n)D’]P.,(k, x) 


(3) 2 —n+1 


+ n| N’ += "| DP,-(k, 2). 


These are not explicitly given in Hildebrandt’s Paper but the method of obtain- 
ing them is outlined there in detail. 

We shall make use of the following lemma which we state without proof. 

Lemma (1). Let P,(x) be a polynomial of degree n. If both P,, and P;, contain a 
factor (x — a)", m <n, then P,, contains the factor (x — a)". 

We also need an expression for P‘?),(k, x). By repeatedly differentiating (2) 
and eliminating P’,(k, x) we get, 


TY 2k —n+7 
Pk, z) = [TL m+1—-% E + a "| P.-11(k, 2), 
(4) “ 
q=1,2,--- (n+ 1), 

Theorem I,. If D is a perfect square, D’ is not a factor of Pasi (k, x), n = 
O£i1.@... 

Proof: Assume D’ to be a factor of P,»;. From (1), D’ is either a factor of 
P,, or of N + (k — n) D’. But D’ is not a factor of N + (k — n) D’ as this 
implies that D’ is a factor of N contrary to hypothesis on (R) that D and N 
have no factorincommon. Thus, D’ is a factor of P, , and by a repetition of the 
reasoning a factor finally of P; , which as it was just pointed out, is impossible. 

Theorem Iz. Set D = (ax + B;)(aox + Be), D not a perfect square. If 
aw + B;,7 = 1 or 2, is a factor of P,, then (ai + B;)* is a factor of Pnsg, 
q=1,2,3,.--- 

Proof: From (1), aw + 8; being a factor of P, and D, is also a factor of 
Psi. From (2), ast + B; is a factor of P),4,. From Lemma (1) it follows 
that (a + £6;)° is a factor of P,4;. Continued repetition of the reasoning 
establishes the theorem. 

Corollary. If both ox + B, and asx + f: are factors of P,, , then D* is a factor 
OC Peiwts 

Theorem I;. Assume D of the same form as in Theorem I,. If ax + 8, 
t= 1 or 2, is a factor of P+, and no higher power of ax + 6; is such a factor then 
at + B; is a factor of N + (k — n)D’. 

Proof: From (1), aw + 6; being a factor of P,.; and of D is also a factor of 
either N + (k — n)D’ orof P,. But aw + 6; a factor of P, requires, from I2 , 
that (a,x + 8,)° be a factor of P,,,, contrary to hypothesis. Thus, aw + 6; isa 
factor of N + (k — n)D’. 
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Corollary. If (ax + Bi)(aex + Be), (a1, a2 ¥ 0), ts a factor of Pri, and no 
higher power of either a,x + Bi or asx + Be is contained in Py41 then N + (k — n) 
D’'=0. For from I; ,N + (k — n)D’ contains (a,x + B1)(aex + Be) as a factor 
which implies N + (k — n)D’, being linear, vanishes identically. 

Theorem Ix. If (ax + Bi)" and no higher power of ax + 8; is a factor of 
Priq— then ax + B; and no higher power of ax + 6B; is a factor of P, . 

Proof: Let us write, 
(A) Paso = (ai + Bi)? Oni, bn1 = a polynomial of degree < n — 1 which 
does not contain the factor aw + 8;. Taking the (q — 1)** derivative of (A) 
by Leibnitz Theorem, we get, 


(B) Pri = > ~ ') & (a:x + B;)! os bn 
n+q-—l = Z dx t a/ dxw-i n—1- 


On setting g = q — 1 in (4) there results, 


q—2 ‘ : 
(C) Pes = T@t+q-1-% [7 ea t-3tit2 "| P,. 
i=0 





2 


From (B) we see that aj + 6; is a factor of P\72,. No higher power of 
a,x + 6;is such a factor. From (C) our theorem now follows. 

Corollary (1). Under the hypotheses of Theorem I,, ax + 8B; is a factor of 
N+ (k—n+1)D’. This follows at once from J, and 73. 

Corollary (2). If D* = (ax + Bi)* (asx + Be)’, (a1, a2 ¥ 0), ts a factor of 
Prsq—-1 and no higher powers of either a,x + Bi or aex + Be are factors, then N + 
(k —n+1)D’ =0. For the linear expression N + (k — n + 1)D’ contains, 
from Corollary (1), the quadratic factor (a;2 + 81) (aor + Be). 

The following lemma can be easily established and is given without proof. 

Lemma (2). Assume D of the same form as in Theorem Iz. Then there is only 
one value of s for which N + sD’ contains a,x + B; as a factor. 

TheoremI;. Assume D of the same form asin TheoremI,. If N + (k — n)D' 
contains ax + B;,7 = 1 or 2, as a factor, then P,+, contains ax + B; and no 
higher power of aix + B; as a factor. 

Proof: From (1) we see that P,.; contains a;z + 8B; at least to the first power 
as a factor. Again from (1), if Pn; contains a higher power of a:x + 8; as a 
factor, this means that both P, and P’, contain a,x + 6; at least to the first 
power as a factor and from Lemma (1) it follows that P, contains a;z + 8; at 
least to the second power as a factor. By corollary (1) from Theorem J, it 
follows that a;z + B;is a factor of N + (k — m)D’ for n, < n, contrary to Lemma 
(2). 

Theorem Ig. If ax + By, and asx + fe are factors of N + (k — 1)D’ and 
N + (k — m)D’ respectively, (a; , a2 ~ 0), then P, =O, np >m+m. 

Proof: From Theorems I; and I2 we see that (az + B:)"* (asx + Bo)”, of 
degree n; + me, is a factor of P,.4n,, of degree m2 + m, at most. Similarly, 
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(at + B:)"?** (aox + Be)"'**, of degree no + m + 2, is a factor of Pryiniis, of 
degree Nm. + m + 1 at most. This implies P,,4n,4. = 0. Hence, P, = 0, 
u>m+m. Infact, (1) shows that P, = 0implies P, = 0,» > u. 

Theorem I;. Assume D of the same form as in Theorem I,. Then Pry, = 9, 
P,, 4 0, implies either N + (k — m)D’ = 0, m < n, or there exist two values of 
m, (m,, M2), such that N + (k — m)D’, N + (k — m2)D’ contain as factors 
art + Bi and aor + Be respectively, (m, , m2 <n). 

Proof: Setting P,i1 = 0 in (1) gives, 


(1°) [VN + (k — n)D'] P, + DP’, = 0. 


If P, = const., 1° shows that N + (k — n)D’ = 0 and our theorem is verified. 
Suppose P,, const. We get from (1°), 
’ [N + (k — n)D']P., 


P, = -—— 


. Thus, 2 is a factor of the numerator, and our theorem now follows from Corolla- 
ries (1) and (2) of Theorem J, . 

TheoremI,. If N + (k — m)D' €0,m = 1,2, --- n, and if N + (k — m)D’ 
contains neither ayx + By , nor aex + Be as factors, then P,,.; and D have no factors 
in common. This follows at once from Theorems J, and J, which constitute a 
necessary and sufficient condition that P, and D have factors in common. 

Theorem I,. If N = const. and if D 1s linear, all P,, are constants, n = 1, 2, 3, 

This follows directly from (2). 

Theorem I. If N’ + = = = 
which are not zeros of D are simple. 

Proof: Suppose P,, has a multiple zero x = @ which is not a zero of D. Then 
(1) shows that a is a zero of P,4;. From (2), a is a zero of Py4,;. From 
Lemma (1), @ is at least a double zero of P,4;. Furthermore, (3) shows that a 
being a double zero of P, and of P,4; is also a double zero of P,,. By a con- 
tinued application of (3), it follows that a is a double zero of P; which is impos- 
sible since P, is of degree < 1. 


D” £ 0, m = 1,2, --- (n — 1), all zeros of P,, 


II. Concerning the Zercs of P,,(k, x) 


The polynomials P,,(k, x) are aefined by Hildebrandt’ as follows: P,(k, 2) = 

—_- 

dx” 
equation 


D*y where y is a non-identically vanishing solution of the differential 


1 dy do + ax _wN 


ydx by +hr+h2 D' 


°L.c. pp. 400-401. 
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The Jacobi Polynomials are defined as follows: 
J Az, a, 8) a (1 oa x)? =. [ate a os. a, 8 


real. It follows that J,(x, a, 8) is a special type of P,(k, x) with N = (—8—a) 
zt+a,D=x(1 — 2x),n = k + 1, whence, 


N’ = —f-a, D=1—27, D” =-2; D@) = D(i) =0, 
P,(k, 2) = N+ kD’ = Ofor 


a 4 k ’ 
L= —— Pi(k, x) = —B — a — 2k. 
- + B 4 Qk? 1\fvy ) B a 
In determining the number and location of the real zeros of the Jacobi Poly- 
nomials we employ the following notations: 


Pk, x) = Oforz = ajx,;, 1=1,2,---k+1;k =0,1,2,---37 =1,2,---74. 
Qi,k,j < Qi.-l,ky7 


2c — 


@ = N’ 
+ 2 


D” = —B—a—2+n, n= 1,2,---k, 


np = (NW + & — aD = a + & — 20), 
paW+0-—- oP. = «+0 —- 6 -— wo. 


We proceed to determine the number of real zeros of the Jacobi Polynomials 
on the intervals (— «, 0), (0, 1), (1, ©) into which the zeros of D divide the 
x axis. The proofs proceed by mathematical induction. We first determine 
the location of the real zeros of P,(k, x), n = 1, 2,---k + 1, by successive 
applications of (1) and (2). We then use the relation P;.4: (k, x) = Js (2, a, B). 

Several cases concerning possible values of a and 8 should be considered. In 
order to bring out the method of procedure only two such cases will be fully 
discussed here. The results for other possible cases will be merely listed. 

Ai: a<0,8 <0,|a| <|8|,a, B, a + B not integers. 
Let k, be the greatest integer contained in a, 

‘cc ko ‘ec ‘cc ‘ec ce éc ‘ec B, 


‘* ks be the greatest integral value of k for which a + 8 + 2k < 0. Then 


0 < ky < ks < ke ° 


4 In the case a, B > 0 these zeros all lie, as is known, on (0, 1). 
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An: 0<kK<k. We then have 0 > 0,4 <0,» >0,0 < ay <1,Pi, > 0. 
Then Ji4i(x, a, B) has a 7 zeros in 0,1. These are the only real zeros. 

Proof: Consider first P\(k, x). Its only zero is at a;,4,, where 0 < ax. <1. 
Furthermore, P} > 0. Also P,; > Oforz > a4. and < Oforz <a 4... From 
(1) we see that P2(k, a1,.,1) > 0, (since Pi(k, a1,%,1) = 0, D(aix.1) > Oand P; > 0). 
From (2) it follows that P2(k, x) < Ofora < a1.%.1, P2(k, a,%1) = 0, Pa(k, x) > 0 
for x > a1,x,. These conclusions follow from remarks concerning the sign of 6, 
the fact that Py(k, o1,x,1:) = 0, and from remarks concerning the sign of P; to the 
left and to the right of x = a,,,. Thus, P2(k, x) > 0 for all real x and hence 
has no real zeros. By employing (2), it is now evident that P3(k,z) >0. From 
(1) and remarks concerning p» and »v we see that P;(k, 0) < 0 and P;(k, 1) > 0. 
Thus P;(k, x) has a single real zero a3,4,0 < a3,z, <1. The reasoning from 
P; to F, is analogous to that from P; to P2. By continuing this procedure we 
finally conclude that Pxi:(k, x), (= Jess (x, a, B), has but one real zero, (in 0, 1), 
if k is even and no real zeros if k is odd. 

Avniki <k < ks. Stk =k +9,q = 1,2,---,ks —k,. Here @> 0, 
up>0,n=1,2,---q—lw<On=¢,qat1,---,aqth. v>O, ana <9, 
Pi(k, z) > 0. Jz, + 9 +1 (2, @, B) has q distinct zeros in (— ©, 0) and 
(1)"* + (-1)"" 

oe 

Proof: First consider the sequence P,(k, x) n = 1, 2, --- g, since the conditions 
on 6, wu, and »v do not change over this range of n. Now Pi(k, a1.) = 0, a4. < 
0. Furthermore since P; > 0 we have P,; > 0 for z > a4, and < 0 forz < 
a. Pass now to P2(k, x). Since D(ay,.1) < 0 and P; (k, a1,x,1) > 0, we see 
from (1) that P2(k, a1) <0. Moreover (2) shows P3 (k, a1,4,1) = 0, P2 (k, 2) 
<0 forz < a, and > Oforz > a x,. Thus Po(k, 2) < 0 and a relative 
minimum at x = a. Since | Po(k, +0)| = ©, we see that P2(k, x) has two 
real zeros of which the left most, a2.,,, isin (—*«,0). Again » > 0 together 
with (1) assures Po(k, 0) > 0. Thus ae,x.2 is in (a1,z,, 0), hence in (— ~, 0). 
By continuing this reasoning on the successive P,(k, x), n = 1, 2, --- q, we 
conclude that P,(k, x) has q zeros in — ©, 0 and P,(k, ax.) < 0. 

Next, consider the sequence P,(k, x), n =q+1,q+2,---qth4+1. 
Over this range of n we have 0 > 0, nu < 0,» > 0. From what has just been 
shown, P,(k, ax.) = 0, — © < agri, < 0,2 = 1,2,---q. Also Pi(k, a,x), 
1 = 1, 2, --- q, is alternately negative and positive. Suppose q odd, (similar 
reasoning holds for qg even). Thus, we suppose P/(k, ag.41) < 0, Pi(k, a¢.k.¢) < 
0, P.(k, x) > Oforz < agi, and < Oforz > agi,g. (1) shows Posilk, a¢,x,:), 
7 = 1, 2, --- q, to be alternately positive and negative. Thus, the zeros a,,x,; 
are separated by gq — 1 zeros of Pyis(k, x). Since from (1), Posi(k, ag.c) > 0 
and from (2) Piii(k, x) > O for x < ayx1, there exists a zero a@oi1,4,1 in (— &, 
ax). Thus far, we have established the existence of q zeros of P,.:(k, x) in 
(— ©,0). q being odd, we have from (1), Paii(k, a@¢,x%,¢) > 0. Also from (2), 


zeros inQ,1. These are the only real zeros. 
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Piii(k, x) < Ofor x > ag,¢. Again from (1) and assumptions regarding yu and 
v it follows that Piii(k, 0) > 0, Paui(k, 1) < 0. Thus, Pyii(k, x) has a zero 
@o+1,k,q41 N (0, 1). There being no extrema for P,.:(k, x) other than the a,,,;, 
a = 1, 2,---4q, (as (2) shows), we have thus proved that P.,:(k, x) has q 
distinct zeros in (— ~,0) and a single zero in (0,1). Reasoning similarly from 
Polk, x) to Pay2(k, x) we establish the existence of q distinct zeros ags2,4,;, 
a= 1, 2,--- q, in (— ©, 0) with agona In (— ©, @oitea) aNd agieks, 7 = 
2, 3, --- g, separating agiix,i, ? = 1, 2,---q. From (1) we see that Py, o(k, 
Qott,kg) < O and Poro(k, agy1,k,.q11) < 0. The only extrema of P 42(k, 2), 
(as (2) shows), are located at agiiz,i, 7 = 1, 2,---q@ + 1. Again, by (2), 
Phas(k, x) < O for x > a¢41,x,941 ; hence there can be no real zeros of P42 except 
the q zeros in (— , 0) already found. The reasoning from P,42 to Pas is 
similar to that from P, to Pgi:. Thus, Pose,41 = Je,4¢41 has q distinct zeros in 
(— ©, 0), together with one zero in (0, 1) for k, even. For k, odd, there are q 


distinct zeros in (— «©, 0) only. The results are the same whether q is odd or 
even. 


The results for the remaining sub-cases under case A, are given in the table 
which follows. For completeness, the results for cases Ay; and Aj: are included 
in the tabulation. A few words of explanation are necessary to clarify the 
conditions under which the various sub-cases in the table occur. Let | a@| = 


ki+4q4,|B8|=kt+h,hq<1. Ifg+h<1,then|a+8| = k, + ke and we 
have either, 


Ays; : ky + ke even, 2k3 = ky + ke = kg — ky = ke — kz 

Aize : ky + ke odd, 2k3 = kp + ke — 1 = kz — ky = ke — keg — 1. 
Again ifl <q+h < 2,then|a+6| = ki + ke + 1 and we have either, 

Aj33 : ky + ke + 1 even, 2k3 = ky + ke + 1 = kg — ky = ke — kg + 1. 

Aisa : ky + ke + 1 odd, 2k3 = ky +*ke = kg — ky = ke — ke 
In cases Ay and Ais; we assume | a + B| = kj + ke + p, p < 1, while in cases 
Aye and Aye ,|a + B| = ki + ke + p,1 < p < 2. The complete results for 
case A, follow. (See page 213.) 

As: a<0,B8 <0,|a| <|8|, a, B not integers, a + 8 = integer. Define k,, 
ke, kj asin Ay. ThenO < kj < ks < ke. In Case Ag , 8 + @ is odd while in 
Case Ag , 8 + aiseven. (See page 214.) 

A;:a < 0,8 < 0,a = —ky, integer, B not an integer,|a|<|B|. Define 
ki, ke, kg asin Ay. ThenO < ki < kj < ke. There are two sub-cases, A3,: the 
greatest integral value of a + 6 is odd, Age : this integral value is even. (See 
page 215.) 

Ay: a < 0, B < 0, @ not an integer, B = —k,, integer,|a|<|B|. Define 
ki, ke, kj asin Ay. Then O < ki < kz < ke. There are two sub-cases, Au : 
the integral part of a + Bis odd, Aw : this integral value iseven. (See page 216). 
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As:a<0,8 <0,|a|<|8|,a= —k, integer, 8B = —k» integer. Define 
ky, ke, kj asin A,;. In cases As; and Age , a + 6 is odd and even respectively. 







Polynomial Range of Sub-Script Zeros in 


—(-»,0) 2=0 £1) 






(1)F + (—0F. 


Asu, Asay Jk413 O<k<h; 0; 0; 9 ; 









Asi2, As2a3 Ski +0413 q = 0,1, 2, --+, ks — ky; q ky +1; 0 












Asis; Jiksta+i3 q 1,2, ---, ke — kg — 1; kg—ki—q; kh+1; 0 




















Aso3; Tkst+q413 q=1,2,--:,ke—kg—1; kg—-hki—q+1; k&&+1; 0 















As Aces; Jkot+a+1 = 0; qg = 0,1, 2, +--+ kt. 








Asis, Aso; Jkitketoti =0; g = 1, 2,3,--- 








If assumptions are identical with those of As; except |a| = ||, then for 
0<k < ky, the results agree with Asu and J;,4¢4. = 0,¢g = 0,1,2,---. 

Ag: a>0,8 <0,|a|>|B8|,B8notaninteger. Let k, be the largest integer 
in B. 





Case Polynomial 


Range of Sub-Script Zeros in 










0,1) (1,0) 
Asi J ie41 0<k<kh 0 sS 
(1)" + (—1)"! 


Ace J ky +041 q= a, 2, 3, oe q 











2 


A; : Same assumptions as in Ag except B = —k, , integer. 








Case Polynomial Range of Sub-Script 


Zeros in 


— 01 r=] (1, 2) 








k k 
—— 0<k<h—-1 0 0 2a 


















An Jir+aut . 0,1, 2, sd q ki+1 0 





As: a > 0,8 < 0,|a|=|8|. J: = a@ and results for J,, n > 1 are 
identical with those in Az and Ag respectively according as £ is or is not an integer. 
Ay: a> 0,8 <0,|a| <|B8|;8,a+ 8B, not integers. 
Let k, be the greatest integer in a + 8B. 
éé ““c“ “cc ‘ec ‘ec “ec 
ke B. 
?_ * o “for which a + B + 2k < 0. 








218 FRANK S. BEALE 


Then 0 < kj < ky < ke. 





Case Polynomial Range of Sub-Script Zeros in 
(—~,0) (,1) (1, ~) 
k+1; 0; 0 
































Ag; = Jks13 O<ck< ks; 


Agai;)  Jkst+qi13 9g = 1,2, ---, kg; ky even; kgs —qt+1; 0; 0 





Ageo;  Jkstqui3 Gg = 1,2, --:, (k3 + 1); kyodd; kj3—q+2; 0; 1 
1)*it¢ —1)*ita 
Ags; Jk +a413.  G = 1, 2, ---, (ke — i); 0; 0; ( tow 
1)*2 + (-1)* 
Ags Jkotant3 «9 = 1, 2,8, +5 0; q; Ors = 





Ay : Same assumptions as in Ag but now|a|=|8|. Then ki = ks = 0, 
J, = a, and results for J, ,n > 1 are the same as in Ag3 and Ag, . 
An : Same assumptions as in Ag except B =-—ke , integer. 





Case Polynomial Range of Sub-Script Zeros in 


(—~,0) (0,1) «=1 (1, ) 





An,1 Same as Agi 
An,2 
An,3 
Aus 





Same as Ago 
















Same as Ag; 











q=1,',3,--:; 0; q3 ke+1; O 


Tica +a413 





Ap: a> 0,B < 0,|a| <|6|, B not an integer. a + B = odd integer. 
Define k, , ke , k3 as in Ag. 
Aj3 : Same assumptions as in A except a + B = even integer. 








Cases Polynomial _____ Range of Sub-Script _ Zeros in — 
(— co ’ 0) 


Ajyo,1, Ais,13 Same as Ag; 


\Snsverti q = 1,2, ---, ks; kk3—qt+1 
Ax,2; 4 ° 
\Jexs+3 = const. > 0; 
(Sesto q=1,2,---,ks +1; k—qt+2 
Ai3,2} } 
| J 2x5 43 = const. > 0; 


Ayw,3, Ais,3; Same as Ags 





Aiw,4, Aj3.4 ; Same as Ags 
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Aw: Same assumptions as in Ay», except B = —ke integer. Cases Aya, 
Awe and Ay3 have the same results as Aj, Aiz2, and Ajs,3 respectively. 
Aus has the same results as Aqy,4. 

Ais : Same assumptions as Aj3 except B = —ke, integer. Cases Ays., Ais, 
and Ais,3 have the same results as Aj3,1, Ais,2, and Aj3,3 respectively. Ajs,4 has 
the same results as Aq4. 

Ai: a = 0,8 < 0, B — not an integer. 

Let k, be the largest integer contained in B. 
‘* ks be the largest integer for which 6B + 2k < 0. 


Case Polynomial Range of Sub-Script Zeros in 
(—~,0) 2w#=0 (0,1) (1, «) 
Awa; Jk+1; O<ck<ks; k; ; ; 0 


ks — q; : : 0; ki even 
ks — 9 +1; ; ; 1; ky odd 


(1)*" + (-1)" 
2 


Ayo; Jkstati1; @ = 1,2, -°:,ki — ks; 
Aws3; Jkitat1; 9g = 1, 2,3, -->; 0; 1: — 


= - k, — odd integer. Define kz as in Aj. 
= —k, — even integer. Define ks as in Aj . 


Polynomial 


Range of Sub-Script Zeros in 


(—«,0) z=0 (0,1) 





Aiz1, Aisi; Same as Aje,1 
Au7,2} JViks+qt}3 1,2,---,ki —ks3—1; k3 — q; 1; 0; 
Ais,2; Tks+a+1} = 1,2,---,k3 +1; ks—qt+1; 1; 0; 


Jik,+1 =0 
Aiz,3, A1s,3; 
Jkyt+qt+1; q = 1, 2,3, +--+; : : - k+l 


Ap: a=0,8 = 0. J; = 0. 
Jixi1 has k — 1 zeros in (0, 1), 1 zero at x =0, 1 zero at x = 1, k = 1, 2, 3, 


From the definition of J,(z, a, 8) it is readily seen that J,(z, a, B) = (—1)” 
J,(1 — x, B, a). Thus, a transformation of x to 1 — z interchanges a and 8. 
The interval (— «, 0) is transformed into (1, ©) and vice-versa. The points 
x = 0 and x = 1 are interchanged. Consequently, in all previous results we 
may interchange properly @ and 8. 

In the foregoing results, the only real multiple zeros that can occur are at 
either z = Oorxz = 1. In the process of determining the degree of multiplicity 
of such zeros use was made of Theorem 2 . 

Points of Inflexion. By taking (4), setting k = n, and replacing N’ and D” 
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by their values for Jacobi polynomials, we get: P4.(n, 2) = (n + 1) (n) 
B+a+n] [8 +a+ n+ 1] Pri(m, x). From definitions of P,(k, x) and 
J (x, a, B) we easily verify that, 


P,n+4q4,2) = J.(z,a+q+1,8 +q +41), whence, 
Jn(x, a, B) = (n+ 1) (n) [Bat nj [Btatnt 1] Jnr (z, a + 2, B+ 2). 


We conclude that if neither a + 8 + n nora + 8 + n+ 1 vanishes, the points 
of inflexion of J4:(x, a, 8) are at the zeros of odd order of Jna(z, a + 2, B + 2). 

The Degree of J,(x, a, 8). In analyzing the results of cases A; to Ajg inclusive, 
it is noted that in some cases the number of real zeros of J, is less than n. The 
question naturally arises whether the degree of J, is n or less, for then we can 
determine the number of its imaginary zeros. The explicit expression of 
J,(x, a, B) is known from which the degree of J, can be found for various a and 
8B. However, the degree of J, can be found from (4). 

Since Jn4i(z, a, B) = Prsi(n, x), let us replace k by n in (4) and at the same 
time replace N’ and D” by their values for Jacobi Polynomials. Thus, we get: 


q—1 
JS21(x, a, B) = I (n+ 1—2)[-—B —a—n—@]Pr_qiiln, 2), 


(5) 


n=0,1,2,---;¢=0,1,---,(n+ 1). 





We may establish the following results. 

C,) If a + B is not an integer, the degree of Jay: (z, a, 8) isn +1, =), 
iS uss; 

In fact, in order for J ‘a to vanish, we see from (5) that either some factor 
— B — a — n — ti vanishes or P,,_,4:(n, x) vanishes identically. We first show 
that the latter is not possible. Now Pi(n, 2) = N + nD’ = (— B — a — 2n) 
zx+a+n# Osince B + ais not an integer. Consequently, if P,(n, x) = 0, 
u>0,u < n+ 1 there will be a first value of u, (u = v), for which P,(n, x) = 0 
but P,.(n, x) # 0. By virtue of Theorem J; this means that either N + 
(n — p)D’ = [— B-—a— 2An— p)|x+a+n— p=), p < », or else there 
exist two values of p, (pi , pe), such that [— B — a — 2(n—pi)Ja+atn—-p 
and [— B — a — 2(n — po)| x + a + n — pp are divisible by z andl — 2 
respectively, pi, pp < v — 1, pi X po. Since, however, a + £ is not an integer 
we see that, [— B — a — 2(n —-p)]x + a+n— p#0,n and p being integers. 
This eliminates the first possibility that P,(n, 7) = 0,4 <n+1. Again, if, 
[—- B—a—2n— p)|x+a+n— pris divisible by z, we have a + n — pir = 
0 oraaninteger. For (a +n — po) — [8 + a+ 2(n — ~r)] x = (a + mn — pr) 
E _ (a+n— pm) + B+n— pr) 

(a + n — pr) 
po = Oor B, an integer. aand Bare therefore both integers contrary to hypoth- 
esis. Thus, in (5), no polynomial P,- ii(k, z) = 0 and Jnii(z, a, B) # 9. 
Replacing g by n + 1 in (5) leads to, 





x | to be divisible by 1 — x requires B + n — 
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ISA? (@, a, 6) = Im +1- 8) [-8-a—n—a Pon, 2), 


(6) 
n=0,1,2,---. 
Thus J$%"” ¥ 0, (since Po(n, x) = 1 and no factor — 8 — a — n — ican vanish) 
and the degree of Jn4; is precisely n + 1. From similar reasoning we prove: 
C.) Ifa+ 6 > Othe degree of Jizsisn +1,n =0,1,2,---. 
C;) If a + B = 0, then (I) J: = a and (IL) Jay: is of degree n + 1, n = 1, 
S Bissx., 
C,) Ifa+ 8 = —M — integer, M > 0, B, a not integers, then, 
(1) For n < M, the degree of Jn4; is min. (n + 1, M — n). 
(II) n = M, Jays = const. 
(III) n > M, the degree of Jnsi isn + 1. 
C;) Ifa + 8 = — M — integer, M > 0, a, 6 integers, a > 0, 8 < 0, then, 
(1) For n < M, the degree of Jn4; is min. (n + 1, M — n). 
(II) n = M, Jays = const. 
(IIT) n > M, the degree of Jinsi isn + 1. 
Cs) If a+ 8 = — M — integer, M > 0, a = —k,-integer, 8 = —k,-integer, 
ky a ke then, 
(I) For n < ke, Jnsi is of degree n + 1. 
(II) n > ke ’ dost = 0. 
C;) Ifa + B = — M — integer, M > 0, a = 8B = —k,-integer, then, 
(I) For n < ky, Jny is of degree n + 1, 
(II) n>k, Janu = 0. 
The Laguerre Polynomials. These are defined as follows: 


L, = L, (a, a) = xe .. [era], n=0,1,2,---; 
dx” 


a — real. We see that L, is a special case of P,(k, x) with N = — 2+ a, 
D=2,n=k+1. It follows that@ = -—lru=atk—n,am=a+hk, 
and Pi(k, z) = 1. These can be used in determining the location of the real 
zeros of L, , as was done for J,. The discussion here is somewhat simplified 
since L, has but one parameter, a, and the z-axis is divided by the zeros of D(z) 
into two segments only, namely, (— ~«, 0) and (0, ~). 

The following results are easily obtained. 
Bi: a > 0, L(x, «) has n distinct zeros in (0, ©), n = 1, 2, 3,---. This 
result is well known. 
Be: a= 0. Ln4yi(x, a) has n distinct zeros in (0, ©) and a simple zero at z = 0, 
n= 0,1,2,---. 
Bs: a < 0, a, not an integer. Let k; be the largest integer contained in a. 


k k 
(I) Lysi(x, aw) has = <2 zeros in(—»,0),0<k<hkh, 
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ky ky 
(II) Li,+ci1(2, &) has q distinct zeros in (0, ©) and eon zeros in 
(— 0,0), 9 = 0, 1, 2, aor 
Bu: a <0,a= —k, — integer. 


i dete ite 2 o. ET tthe eee eh. 


(II) Li,+q4:(@, a) has q ne zeros in (0, ©) and a zero of order k; + 1 at 

z=0,¢=—0,1,2,---. 

The Degree of L,(x, a). We show first that here Pu(n, x) ¥ 0, u = 1, 2,.-. 
n+ 1. By definition, Pi(n,z) = N+ nD!’ = -xt+a+nF0O. Let us 
rewrite (2) for our present situation thus: 

(2°) Pi(n, z) = —uP,-4(n, x). If, now, P,(n, z) = 0, then from (2°) it follows 
that P,1(n, x) = 0. Continuing this reasoning, we finally arrive at a contra- 
diction, namely, P,(n,z) = 0. If in (4) we set g = n + 1 and replace N’ and D” 
by their values we get: 
















Lisi (2, @) = (-1)""(n +1)! Pon, z) = (-1)""(n + 1)! 


Hence, Lny: is of degreen +1. Note that this holds regardless of the value of 
: a contrary to what was found for Jacobi Polynomials. 
Points of Inflexion. By a procedure analogous to that used for Jacobi Poly- 
nomials we can show that the points of inflexion of L,4;(z, a) are located at the 
zeros of odd order of Ly_;(x, a + 2). 
The Polynomials P,(0, x). If we set k = 0 in (1), (2), and (3) we obtain the 
following relationships for P,(0, 2)” = P,(x) = P 


(7) Pnss(z) = [N — nD’) P,(x) + DP, (2). 


(8) Pras(z) = (n+ 1) | W” - * "| P,(2). 








(9) Prii(z) = [N — nD'] P,(2) + n(w N’- 1 D’) DP,,._(2). 


Theorems J, to I) inclusive, with k = 0, hold for P,(x). In addition, the 
following theorems hold for P,, . 


Theorem H,. Suppose N linear and D(x) > 0 for all x. Furthermore, let 
N’ — 5 D” < 0,m = 1, 2,3,---. Then P, has n real, distinct zeros which 














separate the zeros of Pn. 

Proof: Denote the zeros of P, by ani, 7 = 1,2, --- n, anys < @n,iz1. Suppose 
N’ > 0. N being linear has a single zero a. Furthermore, since P; = Ni, 
then P; < 0 for z < ay, and > Oforz > a,. We passnow to P,. From (7), 
we see that P2(a) > 0, (since D > 0 and P; > 0). Also (8) shows P2(x) > 0 





5 E. H. Hildebrandt, loc. cit. pp. 399. 
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forz < anand < O0forz > a,. This follows from what was noted concerning 
m 
. 
Dp” <0. Thus, there exists a zero of Pz in (— ©, ay) and a zero in (ay, ©) 
and our theorem holds form = 1. Assume that the theorem is true for n = h. 
The sequence P;(a,,:), 7 = 1, 2,---h, is alternately positive and negative. 
Since, from (8), the only extrema of P),4; are at a,,;,7 = 1,2, --- h, we conclude 
that there are h — 1 zeros of Pi; separating the a,,;,7 = 1, 2,---h. Since 
Pi(an1) > 0 we conclude that P, < 0 for z < a,,. This fact, combined with 
(8), shows Pyii(xz) > 0 for t < ana. Phsa(ana) being positive, it follows that 
there exists a zero of Pri; in (— ©, a). Similar reasoning establishes the 
existence of a zero of Pi; in (a,,,, ©). Our theorem is thus established for 
N’>0. The case N’ < 0 can be similarly treated. 


Theorem Hz: If D(x) > 0 for all x, D” <0, N’ — 5D" <0,N’ =0,N £0, 


the sign of P, for x > ay and x < ay, , together with the hypothesis that N’ — 


then P, , n = 2, 3, --- , has n — 1 real, distinct zeros which ave separated by the 
zeros Of Pn. 

Proof: Since P; = N = const., we see from (7) that P: is linear. The reason- 
ing of Theorem H, applies where we now start with P2. 

Theorem H3: If D(x) > 0 for all x, except x = B, where D has a double zero and 


if N’ # 0, N’ — 5 D” <0,n = 1, 2, 3, --- , then P, has n real, distinct zeros 


which separate those of Pn. 

Proof: Theorem J, with k = 0 assures us that P, and D have no zeros in 
common. The proof now follows the line of reasoning of Theorem H, . 

Theorem H,: If D(x) > 0 for all x except x = B where D has a double zero and 


if N’ = 0, N #0,N’ — 5D" < 0, m = 1, 2,3, --- , then P, has n — 1 real, 


distinct zeros which separate those of Pasi ,n = 1, 2,3, ---. This theorem follows 
from H; as did Hz from H,. 


Points of Inflexion. Setting k = 0 in (4) leads to, 
Py. = (n +1) (n) Ei “ "| Ei ~s 5 : D” | Pia 


This shows, under the assumptions of Theorems H, to H, inclusive, that the 
points of inflexion of P,,; are at the zeros of Pp_;. 

Hermite Polynomials. Theorem H, and statement immediately above con- 
cerning points of inflexion apply directly to Hermite Polynomials where N = —z 
and D = o’. 
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THE SIMULTANEOUS COMPUTATION OF GROUPS OF REGRESSION 
EQUATIONS AND ASSOCIATED MULTIPLE CORRELATION 
COEFFICIENTS 


By Paut S. Dwyer 


1. Introduction. The need sometimes arises for the prediction of a number of 
different variables from a given group of so-called fundamental variables. In 
the work of college prediction, for example, one might desire regression equations 
predicting certain measures of college achievement (e.g., first semester average, 
first semester English grade, first semester mathematics grade, number of hours 
of A received during first semester, etc.) on the basis of a number of other factors 
(e.g., high school record, score on American Council on Education Psychological 
Examination, score on some standard English achievement test, score on some 
standard mathematics achievement test, etc.). It is the purpose of this paper 
to show how the regression coefficients and the associated multiple correlation 
coefficients can be obtained simultaneously. The essence of the method is a 
simple device by which one solution of general normal equations may be made to 
serve for all cases. 


2. The normal equations. Let 2, x2, x3,---2,, be the so-called funda- 
mental variables and let x, be the predicted variable. The normal equations 
are computed by standard methods which result in one of the three types. 

Type I. Normal equations for determining bp , b; , be , bs , - 


bon aa bi D2 oe be=ze + b3D23 + 
1 2 . ~ 

bola, + bi 221 + betaixe + b32aix3 + 
5 2 . 

bo2Xe + by D222 + bel xe a bs Xe23 oS 


Type II. Normal equations for determining b; , be, 
Xi = Li — M,, 


b2Zi 0 + beLEF2 + bs2HZ,%3 + 


ie v2 
bi DFoT, + DoDZXo 
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Type III. Normal equations for determining 6; , B2 , Bs, --- , Bn. 
Bi + Tbe + 11383 + 
Tabi + Be + 12383 + 


The three types are special cases of the general 
dyyr + diye + disy3 + 
doyi + dooye + desys + 
dsiyi + dsey2 + dssys + 


where y; are the regression coefficients and d;; = dj; . 

The methods described in this paper are applicable to the general case 
hence to each of the three particular types. 

In examining the normal equations, it is noticed that the first n terms of each 
equation are completely determined by the n fundamental variables. The 
equations, aside from the last terms, are identical no matter what variable is 


predicted. It is only necessary to devise a technique for separating the con- 
tributions of the dj, terms. 


3. Solution by determinants. One method utilizes determinants. The 
value y; is expressed in terms of a determinant involving a column with entries 
dy, dex, dsz, --- , dnx. The determinant is expanded in terms of this column. 

Specifically, let D be the determinant of the coefficients of the y; and let D;; 
be the cofactor of any element d;; of D. Then 


D = > Duds 
i=1 


= = (Dudy + Deda, + Dads, +....+ Dard +....+ Dri dne.) 


= D (Dio dy, + Doe dy. + Ds2 dg. +....+ Djedj, +.... 


= 5) (Din ds + Don dx + Dsndse +...-+ Dindix +....+ Dan dne-) 
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; Di: ‘ ‘ 
It is only necessary to compute - D to find the coefficient of d;, in the expansion 









of Yi. 
An illustration is given. The normal equations are 


Bi + .3300 B. + .2100 63 — ry = 
.3300 6; + Be — .4800 B3 — ro, 
.2100 Bi — .4800 Bo + Bs — 13% 















from which at once 
a * (.7696 ry. — 4308 rx, — .368473,) 


1 


— D (—.4308 re + .9559 re, + .5493 rei) 


— : (—.3684 ru + .5493 rx + .8911 73.) 


and also 
D = .550072 = (1.00)(.7696) + (.33)(—.4308) + (.21)(—.3684) 

= (.33)(—.4308) + ( 1.00)(.9559) + (—.48)(.5493) 

= (.21)(—.3684) + (—.48)(.5493) + ( 1.00)(.8911) 


so that 















1.3991 ri. — 





Bi = 
Be 
Bs = —.6697 Tik + .9986 Tok ob 1.6200 T3k - 


.7832 1) ie .6697 3k - 





— .1832 T1k + 1.7378 Tok oa .9986 13k - 





It is only necessary to insert any given values 7% , 72: , 73; , to obtain the coeffi- 
cients of any specific regression equation. 


4. Solutions without determinants. Theoretically the solution by deter- 
minants is excellent but as the number of variables increases the work of com- 
: n(n + 1) ,. 
puting the n’ cofactors | the = s different cofactors becomes enormous. 
We desire a technique for separating the contributions of the last terms when 
determinants are not used. This can be accomplished by using a separate 
column for each dy. Before algebraic manipulation, the value dj, is factored 
from the column and, after manipulative solution is complete, the multiplication 
by dy. is carried out. 
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As an example consider the normal equations 
Bi + ribo — TR = O 
Ba + Be — rx = 0 
where 712 = 72 = .3300. Then the normal equations may be represented by 


rows (1) and (2) of Table I. 


TABLE I 


Tik 


Row Operation | B: | Bo 


(1) 1.0000 .3300| — 1.0000 
.3300' 1.0000, 

— .3300 times (2) | — .1089) — .3300) | ‘ 

(1) + (3) | .8911) | —1.0000, —.3300 

| — (4) divided by .8911 | —1.0000 | 1.1222) — .3703 

| — .3300 times (5) | .3300 | — .3703| . 1222 

| — (2) + (6) | | —1.0000, — .3703| 1.1222 


7) 


The four decimal place solution, whose steps are indicated by (3) (4) (5) (6)(7), 
is from (5) and (7) 


Bi 1.1222 nk — 3103 Tok 
Be = —.3703 ry, + 1.1222 re 


This device may be combined with most of the standard methods of solving 
normal equations. 


5. Combination with Doolittle method. Especially to be recommended is a 
combination of this device with the Doolittle method which is recognized as a 
most efficient method of solving normal equations in from five to ten variables 
[1] [2]. One of the advantages of the Doolittle method is that related multiple 
regression coefficients may be obtained from the same forward solution, though 
additional back solutions are necessary [3]. 

The problem which led to the development of this technique was the simul- 
taneous prediction of scores on various occupations covered by the Strong 
Vocational Interest Blank from the scores on a few fundamental occupations. 
A multiple factor analysis revealed that five basic factors account for most of the 
scores. Five occupational scores, serving as approximations to the five basic 
factors, were used as the fundamental variables and the other scores were 
predicted from them. 

As an illustration of this prediction technique combined with the Doolittle 
method, I have selected three test scores as fundamental since the solution based 
on them shows all the steps of the Doolittle method and is shorter than the five 
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variable problem. 
for problems involving three variables. 
presented in Table II. 
determinant solution of section 3. 














algebraic operation. 











Tik , Tex , T3x kept in separate columns. 















































































































































Bi = 
Be 
Bs 


1.3990 ry. — 

















— .6697 ru. + 


TABLE II 


Generalized Doolittle Presentation 


DWYER 


Actually, solution by determinants (section 3) is advised 








The steps of the Doolittle solution are 


The results should be compared with those of the 


The first column indicates the row and the second the description of the 
The next three columns are the standard columns of a 
Doolittle presentation with the conventional elimination of the lower left entries, 
The next three columns carry through the Doolittle method with the values 
The last column is an adaptation of the 
conventional summary check column of the Doolittle solution. 





Row Operation Bi B2 Bs Tik Tok T3k | S 
(1) | 1.0000) 3300, .2100 — 1.0000 | | .5400 
(2) | 3300} 1.0000} — .4800 —1.0000) | —.1500 
(3) | 2100, —.4800, 1.0000 ‘1.0000, — .2700 
(4) | Repeat (1) 1.0000 .3300, + .2100|—1.0000) | 5400 
(5) | Negative of (4)  |—1.0000) —.3300, —.2100| 1.0000 — .5400 
(6) | Repeat (2) | 1.0000, —.4800) |—1.0000) ~ 1500 
(7) | —.3300 times (4) | —.1089 —.0693| 3300 | — .1782 
(8) | (6) + (7) | 8911 —.5493) .3300 —1.0000| — 3282 
(9) | —(8) divided by 1.0000 .6164) —.3703, 1.1222! 3683 
| 9911 | | | 
emanieigreseu or setiennsitnieetrenh ection feminism tn pens resentments 
(10) | Repeat (3) | | 1.0000 | 1.0000 — .2700 
(11) | —.2100 times (4) —.0441) 2100) | | —.1134 
(12) | .6164 times (8) | | | —.3386) 2034) — .6164| | — .2023 
(13) | (10) + (11) + (12) | | | .6173)  .4134) — .6164/—1.0000, — .5857 
(14) | —(18) divided by | '—1.0000, —.6697| 9985) 1.6200] 9488 
| .6173 | | | | 
(15) | .6164 times (14) | | | —.6164) —.4128, .6155) 9986, 5848 
(16) | (9) + (15) | |—1.0000 | —.7831| 1.7377; .9986| 9531 
(17) | —.2100 times (14) 2100, 1406) — 2097) —.3402, —.1992 
(18) | —.3000 times (16) | | 3300) .2584| — .5734| —.3295| —.3145 
(19) | (5) + (17) + (18) |—1.0600) | 











1.3990| — .7831| — .6697 —1.0537 


The general solution is read from rows (19) (16) (14) and is 
7831 To — 
— .7831 ik ote 1.7377 Tok + 


.6697 13k « 
.9986 T'3k « 
9985 rox + 1.6200 rz; . 
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which agrees, aside from the last place, with the result of the solution by de- 
terminants. 


It is wise to check in the original equations (1), (2), (3) as soon as any 8; is 
found. Row (14), for example, should be checked by showing 


(—.6697)(1.00) + (.9985)( .33) + (1.6200)( .21) = .0000 
(—.6697)( .33) + (.9985)( 1.00) + (1.6200)(—.48) = —.0001 


(—.6697)( .21) + (.9985)(—.48) + (1.6200)( 1.00) 1.0001 
The same should be done with row (16) as soon as it is computed. Row (19) 
should be treated similarly. . 


6. Many regression equations. If large numbers of regression equations are 
to be generated (the Strong Vocational Interest Study had 29 dependent va- 
riables), the following technique is suggested. Make a table with columns 
Te, To , ete. and use the rows to indicate the different values of k. On another 
slip of paper insert the general values 6; , B2, Bs, --- Bx in successive rows so 
that a folding of the paper will bring any general 6 expansion in conjunction 
with the r’s of any test, k. The scheme is illustrated in Table ITI. 


TABLE III 


Occupation Tik | Tok 3k 


| Teacher 1.00 33) .21 
Physicist 33 1.00; —.48| 
| Office Worker .21; —.48 1.00) 
Doctor AG 79, = —.52 
Lawyer | —.02) .165 —.59| 
| Engineer 16 .78 — .02; 


1.3990) —.7831) — .6697 
| 








—.7831| 1.7377|  .9986 


| | 
| —.6697, 9985, 1.6200 


10 | Mathematician .46 .96 — .49 
| etc. 





Thus, for the occupation of Engineer, 
B, = 1.3990 (.16) + (—.7831)(.78) + (—.6697)(—.02) 
B2 = —.7831 (.16) + ( 1.7377)(.78) + ( .9996)(—.02) = 1.21 
8; = —.6697 (.16) + ( .9985)(.78) + ( 1.6200)(—.02) = 64 
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The value of the multiple correlation coefficient is then computed from the 
formula 


Tk .12—-—n = V Bux Vik + Box. Yer +. see + BnkT nk 
In the illustration above 
rea = W/(—.37)(.16) + (1.21)(.78) + (.64)(—.02) 
.93 


7. Regression equations by deletion. The method of getting related regres- 
sion coefficients and correlation coefficients, described by Kurtz [3], is also 
applicable. Again, a problem involving more than three variables is needed to 
show the real value of the scheme but the technique may be illustrated in the 
three variable case. We wish to find, from the forward solution of Table II, 
the regression equation and the multiple correlation coefficient when the first two 
fundamental variables only are used. We delete all columns involving test 3 
and complete the back solution as indicated in Table IV, which may be viewed 
as a substitute for the last ten rows of Table II. 


TABLE IV 
(See Table IT) 








Row | Operation Bi Be B; Tik T2k T'3k Ss 
(20) | Repeat (9) —1.0000 — .3703| 1.1222 
(21) — .3300 times (20) 3300 .1222) — .3703! | 

| 

| 


| (5) + (21) — 1.0000! 


1.1222) — .3703) 








The results are 
Bi 1.1222 ry. —.3703 ro . 


Be = — .3703 Tik + 1.1222 Tok . 





and these agree with the results of section 4. 


8. The simplified back solution. In every case in which the 6’s have been 
given in terms of r’s the matrix of the coefficients is symmetric (sections 3, 4, 5, 7). 
One wonders if this symmetry is’ generally true and if it holds for normal equa- 
tions of Type I or Type II. 

Determinants are much more useful in establishing general properties, such 
as the one under discussion, than they are in computing the values of regression 
coefficients in the case of a problem involving many variables. We return to the 
determinant notation of section 3. 

In each of the three types, and hence in the general case d;; = d;;so that Disa 
Di; Dx 


D D Hence the matrix of the 


symmetric determinant, D;; = Dj; and 


coefficients of the solution is symmetric. 
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This result may be used (1) to check the expanded results or (2) to eliminate 
some of the work of the back solution. The n coefficients must be recorded for 
8, after which the column indicated by r,, may be dropped. The first n — 1 
coefficients must be computed for 8,1 after which the column indicated by 
r,1,.. May be dropped, ete. The italicized entries in Table II are the ones 
which are eliminated in this way. The remaining coefficients are sufficient to 
completely determine the symmetric matrix. 

The summary right hand check column can not be readily used in the simpli- 
fied back solution but it is hardly to be recommended anyway. Kurtz [3] 
argues against it on the ground that it is not necessary. The essential check is 
to see that each 8 solution satisfies all of the original equations. 


9. Conclusion. This paper provides a technique for the computation of 
general regression equations and shows how the technique may be combined 
with the Doolittle method in providing a practical means of mass prediction. 
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CONSTITUTION 
ARTICLE [| 
NAME AND PURPOSE 


1. This organization shall be known as the Institute of Mathematical Sta- 
tisvics. 
. 2. Its object shall be to promote the interests of mathematical statistics. 


ARTICLE II 
MEMBERSHIP 


1. The membership of the Institute shall consist of Members, Fellows, 
Honorary Members, and Sustaining Members. 
2. Fellows shall be the only voting members of the Institute. 


ARTICLE III 


OFFICERS, BOARD OF DIRECTORS, COMMITTEE ON MEMBERSHIP, 
AND COMMITTEE ON PUBLICATIONS 


1. The Officers of the Institute shall be a President, two Vice-Presidents, 
and a Secretary-Treasurer, elected for a term of one year by a majority ballot 
at the annual meeting of the Institute. Voting may be in person or by mail. 

(a) Exception. The first group of Officers shall be elected by a majority 
vote of the individuals present at the organization meeting, and shall serve until 
December 31, 1936. 

2. The Board of Directors of the Institute shall consist of the Officers and 
the previous President. 

3. The Institute shall have a Committee on Membership composed of three 
Fellows. A+ their first meeting subsequent to the adoption of this Constitution, 
the Board of Directors shall elect three members as Fellows to serve as the 
Committee on Membership, one member of the Committee for a term of one 
year, another for a term of two years, and another for a term of three years. 
Thereafter the Board of Directors shall elect from among the Fellows one 
member annually at their first meeting after their tlection for a term of three 
years. The president shall designate one of the Vice-Presidents as Chairman 
of this Committee. 

4. The Institute shall have a Committee on Publications composed of three 
Members or Fellows elected by the Board of Directors. The President shall 
designate a Vice-President as Ex Officio Chairman of this Committee. 
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ARTICLE IV 
MEETINGS 


1. A meeting for the presentation and discussion of papers, for the election of 
Officers, and for the transaction of other business of the Institute shall be held 
annually at such time as the Board of Directors may designate. Additional 
meetings may be called from time to time by the Board of Directors and shall be 
called at any time by the President upon written request from ten Fellows. 
Notice of the time and place of meeting shall be given to the membership by the 
Secretary-Treasurer at least thirty days prior to the date set for the meeting. 
All meetings except executive sessions shall be open to the public. Only 
papers accepted by a Program Committee appointed by the President may be 
presented to the Institute. 

2. The Board of Directors shall hold a meeting immediately after their 
election and again immediately before the expiration of their term. Other 
meetings of the Board may be held from time to time at the call of the President 
or any two members of the Board. Notice of each meeting of the Board, other 
than the two regular meetings, together with a statement of the business to be 
brought before the meeting, must be given to the members of the Board by the 
Secretary-Treasurer at least five days prior to the date set therefor. Should 
other business be passed upon, any member of the Board shall have the right to 
reopen the question at the next meeting. 

3. The Committee on Membership shall hold a meeting immediately after the 
annual meeting of the Institute. Further meetings of the Committee may be 
held from time to time at the call of the Chairman or any member of the Com- 
mittee provided notice of such call and the purpose of the meeting is given to 
the members of the Committee by the Secretary-Treasurer at least five days 
before the date set therefor. Should other business be passed upon, any 
member of the Committee shall have the right to reopen the question at the 
next meeting. 

4. At a regularly convened meeting of the Board of Directors, three members 
shall constitute a quorum. At a regularly convened meeting of the Committee 
on Membership, two members shall constitute a quorum. 


ARTICLE V 
PUBLICATIONS 


1. In the beginning, the “Annals of Mathematical Statistics” shall serve as 
the official journal for the Institute. Other publications may be originated 
by the Board of Directors as occasion arises. 


ARTICLE VI 
EXPULSION OR SUSPENSION 


1. Except for non-payment of dues, no one shall be expelled or suspended 
except by action of the Board of Directors with not more than one negative vote. 
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ARTICLE VII 


AMENDMENTS 


1. This constitution may be amended by an affirmative two-thirds vote at 
any regularly convened meeting of the Institute provided notice of such proposed 
amendment shall have been sent to each Fellow by the Secretary-Treasurer at 
least thirty days before the date of the meeting at which the proposal is to be 
acted upon. Voting may be in person or by mail. 


BY-LAWS 


ARTICLE I 


DUTIES OF THE OFFICERS, BOARD OF DIRECTORS, COMMITTEE 
ON MEMBERSHIP, AND COMMITTEE ON PUBLICATIONS 


1. The President, or in his absence, one of the Vice-Presidents, or in the 
absence of the President and both Vice-Presidents, a Fellow selected by vote 
of the Fellows present, shall preside at the meetings of the Institute and of the 
Board of Directors. At meetings of the Institute, the presiding officer shall 
vote only in the case of a tie, but at meetings of the Board of Directors he may 
vote in all cases. At least three months before the date of the annual meeting, 
the President shall appoint a Nominating Committee of three members. It 
shall be the duty of the Nominating Committee to make nominations for 
Officers to be elected at the annual meeting and the Secretary-Treasurer shall 
notify all Fellows at least thirty days before the annual meeting. Additional 
nominations may be submitted in writing, if signed by at least ten Fellows of 
the Institute, up to the time of the meeting. 

2. The Secretary-Treasurer shall keep a full and accurate record of the 
proceedings at the meetings of the Institute and of the Board of Directors, 
send out calls for said meetings and, with the approval of the President and the 
Board, carry on the correspondence of the Institute. Subject to the direction 
of the Board, he shall have charge of the archives and other tangible and 
intangible property of the Institute. He shall send out calls for annual dues and 
acknowledge receipt of same; pay all bills approved by the President for expendi- 
tures authorized by the Board or the Institute; keep a detailed account of all 
receipts and expenditures, prepare a financial statement at the end of each year 
and present an abstract of the same at the annual meeting of the Institute after 
it has been audited by a Member or Fellow of the Institute appointed by the 
President as Auditor. The Auditor shall report to the President. 

3. The Board of Directors shall have charge of the funds and of the affairs 
of the Institute, with the exception of those affairs specifically assigned to the 
President or to the Committee on Membership. The Board shall have au- ’ 
thority to fill all vacancies ad interim, occurring among the Officers, Board of 
Directors, or in any of the Committees. The Board may appoint such other 
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committees as may be required from time to time to carry on the affairs of the 
Institute. 

4. The Committee on Membership shall prepare and make available through 
the Secretary-Treasurer an announcement indicating the qualifications requisite 
for the different grades of membership. 

5. The Committee on Publications, under the general supervision of the 
Board of Directors, shall have charge of all matters connected with the publica- 
tions of the Institute, and of all books, pamphlets, manuscripts and other 
literary or scientific material collected by the Institute. Once a year this 
Committee shall cause to be printed in the Official Journal the Constitution 
and By-Laws and a classified list of all the Members and Fellows of the Institute. 


ARTICLE II 
DUES 


1. Members shall pay five dollars at the time of admission to membership 
and shall receive the full current volume of the Official Journal. Thereafter, 
Members shall pay five dollars annual dues. The annual dues of Fellows shall 


‘ be five dollars. The annual dues of Sustaining Members shall be fifty dollars. 


Honorary Members shall be exempt from all dues. 

2. Annual dues shall be payable on the first day of January of each year. 

3. The annual dues of a Fellow or Member include a subscription to the 
Official Journal. The annual dues of a Sustaining Member include two sub- 
scriptions to the Official Journal. 

4. It shall be the duty of the Secretary-Treasurer to notify by mail anyone 
whose dues may be six months in arrears, and to accompany such notice by a 
copy of this Article. If such person fail to pay such dues within three months 
from the date of mailing such notice, the Secretary-Treasurer shall report the 
delinquent one to the Board of Directors, by whom the person’s name may be 
stricken from the rolls and all privileges of membership withdrawn. Such 
person may, however, be re-instated by *he Board of Directors upon payment 
of the arrears of dues. 


ARTICLE III 
SALARIES 


1. The Institute shall not pay a salary to any Officer, Director, or member of 
any committee. 


ARTICLE IV 
AMENDMENTS 


1. These By-Laws may be amended in the same manner as the Constitution 
or by a majority vote at any regularly convened meeting of the Institute, if the 
proposed amendment has been previously approved by the Board of Directors. 
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